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Abstract 

In Large-Eddy simulation of particle-laden flow, the effect of the unresolved scales on the 
particles needs to be modelled. In this work we analyse three very promising models, namely 
the approximate deconvolution method (ADM) which was proposed for particle-laden flow 
independently by Kuerten (Phys. Fluids 18, 2006) and Shotorban and Mashayek (Phys. 
Fluids 17, 2005) and two stochastic models, proposed by Shotorban and Mashayek (J. Turbul. 
7, 2006) and Simonin et al. (Appl. Sci. Res. 51, 1993). We present results from a priori and 
a posteriori analysis of these models in isotropic turbulence at Re\ = 52. This data allows 
for a direct quantitative comparison of the models. The analysis shows that ADM always 
leads to improved statistics but that even for high Stokes numbers, the rate of dispersion is 
not predicted correctly by ADM. Concerning the stochastic models, we found that with the 
correct choice of model parameters, the models perform well at small Stokes numbers. On 
the other hand, at high Stokes numbers the stochastic models show significant errors such 
that it may be recommendable to neglect the small scale effects instead of using one of the 
stochastic models. 

Keywords: Large-Eddy Simulation, particle-laden flow, SGS effects, Approximate 
Deconvolution, Langevin model 



1. Introduction 



Large-Eddy Simulation (LES) has become an important tool for the simulation of tur- 
bulent flow. State of the art methods provide reliable results and are capable to tackle 
application relevant challenges. One crucial component for LES is the correct choice of a 
turbulence model, i.e., a model for the effect of the unresolved subgrid scales (SGS) on the 
resolved scales. Such models are herein referred to as fluid-LES models. 

For Large-Eddy Simulation of particle laden flow, an additional model for the effect of 
the unresolved scales on the particles is needed, referred to as particle-LES models. The 



works of 


Yamamoto et al. 


(2001) 


; Armenio et al. 


(1999) 


and Simonin 


(2006 


); 


Yang et al. 


2008 


); 


Marchioli et al. 



(2008) show that neglection of small 



scale effects is not an option. 

Most particle-LES models were developed in a Eulerian-Lagrangian framework, i.e., the 
carrier fluid flow is computed by solving the Navier-Stokes equations and the particles are 
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computed by tracing single particles through the domain. Then, modelling reduces to re- 
construction of small scale effect on a single particle. 

On this basis, a large number of models was proposed. Among these are for example 



the models of S 


imonin et al. 


(1993 


); 


Wang and Squires ( 


1996 


); 


Shotorban and Mashayek 


(2006 


); 


Kuerten 


(20061 


>); 


Amiri et al. 


(2006 


); 


Gobert et al 


(2007 


); 


Shotorban et al. 


;2007); 


Bini and Jones 


(2007, 


2008 


); 


Pozorski and Minier 


(1998) 


Pozorski and Apte 


(200£ 


)), just 



to mention a few. Most of these models are stochastic models, often obtained by extending 
models which were originally developed for inertia free particles in the context of Reynolds 
Averaged Navier-Stokes (RANS) simulations, such as the generalised Langevin model by 



Haworth and Pope (1986). 

A deterministic alternative is the approximate deconvolution method (ADM) for particle 
laden flows (see Kuerten, 2006b Shotorban et al. 2007 Shotorban and Mashayek, 2005). 



ADM is based on an approximate inversion of the LES filter and was originally developed 
in a Eulerian context. 

The present study focusses on three very promising particle-LES models, namely the 



stochastic models proposed by Shotorban and Mashayek (2006) and Simonin et al. (1993) 



and ADM as proposed by Kuerten (2006b) and Shotorban and Mashayek (2005). For all 
three models, the respective authors present some results on the accuracy of their models. 
Their findings are summarised as follows. 



Kuerten (2006b) analysed ADM in particle-laden turbulent channel flow at a Reynolds 
number based on friction velocity of Re T = 150. He conducted an a posteriori analysis for 
particles with Stokes numbers of St = 1, 5 and 25, based on the viscous time scale. His 
results show that ADM significantly improves rms values of the wall normal component of 
the particle velocity. The improvement is greater for high Stokes number than for low Stokes 



number. In addition, Shotorban and Mashayek (2005) found that in a turbulent shear layer 



ADM improves particle dispersion. 



Shotorban and Mashayek (2005 ) analysed their Langevin-based model in decaying isotropic 
turbulence and found that for small Stokes numbers (St < 2.5 based on the Kolmogorov 
time scale at initialisation) the model leads to correct particle dispersion whereas at higher 
Stokes number significant deviations can be observed. 



Fede et al. (2006) analysed the model of Simonin et al. (1993) in forced isotropic turbu- 



lence and found that the model leads to correct kinetic energy for the particles. However, 
their simulations are restricted to St < 5, based on the Kolmogorov time scale. The present 
study shows that at higher Stokes numbers the model does not perform very well. 

Concluding, all published results were obtained on different configurations and are there- 
fore not comparable. In particular, for the Langevin-based models only data at small Stokes 
numbers is published. For all models, the available data density over the Stokes number 
range is not satisfactory. Data rather correspond to probes at specific Stokes numbers but 
from this data no Stokes number dependent behaviour of the models can be deduced. 

The present study aims at a clarification of that issue by providing data which allows a 
direct comparison of these three particle-LES models on a broad range of Stokes numbers. 
The data density on the Stokes number range is sufficiently high to allow the deduction of 
a Stokes number dependence. The testcase is isotropic turbulence at Re\ = 52. All three 
models were originally developed such that they should perform well in that testcase but we 
will show that even by tuning the model constants, the models do not always perform well. 
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Actually in some cases better results are obtained by neglecting SGS effects than using one 
of the stochastic models. 

This paper is organised as follows. Sections [2] and [3] contain a description of the nu- 
merical methods used to compute flow and particle dynamics. Statistics of the single phase 
simulations are also presented in section [2} In section [4] the three particle-LES models under 
consideration are presented and section [5] contains results of an a priori and an a posteriori 
analysis of the models. 



2. Numerical Simulation of the carrier flow 



In the present work we analyse particle dynamics in forced isotropic turbulence by DNS 
and LES. For the simulation of the carrier fluid, we use a second order Finite- Volume method 



together with a third order Runge-Kutta scheme proposed by Williamson ( 1980 ) for advance- 



ment in time. The conservation of mass is satisfied by solving the Poisson equation for the 



pressure using an iterative solver proposed by Stone (1968). More details on the flow solver 



can be found in Manhart (2004). 



The flow is driven using a slightly modified version of the deterministic forcing scheme 



proposed by Sullivan et al. (1994). Sullivan et al. propose a forcing scheme where the energy 
in the spectral modes below a certain wave number K\ is held constant. We additionally 
imposed a lower bound for the forced wavenumbers, i.e., only the modes in a given range 
[fto, K i] are forced. The Reynolds number in our simulations is always Re\ = 52, based on 
the transverse Taylor microscale A and the rms value of one (arbitrary) component of the 
fluctuations u rms . 

In all computations the flow was solved in a cube on a staggered Cartesian equidistant 
grid. The size of the computational box and the cell width was chosen such that all scales 
are resolved, based on the criteria stated by Pope (2000), cf. table IT] 



Table 1: Simulation parameters and Eulerian statistics from DNS of forced isotropic turbulence. 



DNS 

~Re~ x 52 

Number of grid points N 256 3 

range of forced wavenumbers [k , Ki] [0.514, 1.54] /A 

integral length scale Lf 2. 00 A 

time scale of energy containing eddies kf/e 5.15X/u rms 

Kolmogorov length scale t)k 0.070A 

Kolmogorov time scale tk 0.248A/w rms 

length of computational box L 11.9Lj 

cell width Ax 1.3At}k 

filter width A 7 Ax 

kinetic energy of the filtered field kf 0.87 /kf 



The particle-LES models were assessed by a priori and a posteriori analysis. For the a 
priori analysis, we filtered the DNS field u by a box filter Q with filter width A = 7Ax, Ax 
being the DNS cell width. The filtered DNS was sampled on a correspondingly coarse grid, 
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resulting in a field Qu which is comparable to an LES field. The kinetic energy of the filtered 
field k f = (Guf) /2 is 87% of the energy of the unfiltered field k f = (uf) /2, cf. table [l) (■) 
denotes spatial and temporal averaging. 

In the a priori analysis, Qu was used as input for the particle-LES models. Then, the 
models were assessed with respect to the difference in statistics obtained from unfiltered 
DNS and filtered DNS with particle-LES model. 

For the a posteriori analysis, Qu was computed by LES. As fluid-LES model we used the 



Lagrangian dynamic Smagorinsky model proposed by Meneveau et al. (1996). The forcing 
parameters for LES were chosen in the same way as for DNS, i.e., the energy contained in 
the range [k ,Ki] is equal in LES and DNS. Beyond K\, the energy in LES is lower than in 
DNS due to the different grids and the fluid-LES model. With our choice of the grid, the 
kinetic energy resolved by LES kf is approximately equal to the kinetic energy of the filtered 
DNS field kf, cf. table [2J Instantaneous energy spectra E(k) from DNS and LES are plotted 
in figure [l] In addition, a model spectrum proposed by Pope (2000) is shown. All data were 
made dimension free by normalising with DNS quantities. 

Table 2: Parameters for LES of forced isotropic turbulence. 

LES 

Re x 52 

Number of grid points A" 42 3 
cell width Ax 0.567A 
time scale of energy containing eddies kf/e 10A0X/u rms 
resolved kinetic energy kf 0.86kf 




3. Discrete particle simulation 

In this study we consider dilute suspensions of small particles. Thus, effects of the 
particles on the fluid and particle-particle interactions are neglected (one way coupling). 

The density of the particles was set to p p = 1800p where p is the density of the fluid. 
In each simulation the particles were divided in 24 fractions with different diameter d. The 
maximum diameter equals the Kolmogorov length scale. Consequently, the particles can be 
treated as point particles. 
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The particle relaxation time 



p P d 2 



v p l%v 

ranges from t p = O.Itk to r p = 100r^. Corresponding Stokes numbers St 
the Kolmogorov time scale tk range from St — 0.1 to St = 100 . 



— based on 



Based on the works of Armenio and Fiorotto (2001) and Kubik and Kleiser (2004), we 



assumed that in the given configurations the acceleration of a particle is given by Stokes 
drag only, 



dv 

~dt 



cpRe p 
' 24tv, 



u 



(2) 



Here, v(t) denotes the particle velocity and u/® p the fluid velocity at the particle position. 
The particle Reynolds number Re p is based on particle diameter and particle slip velocity 
||u/@ p — v || which leads to a nonlinear term for the Stokes drag. The drag coefficient cd was 
computed in dependence of Re p according to the scheme proposed by Clift et al. (1978). 

The fluid velocity Uf &p must be evaluated at the particle position x p (rj), i.e. Uy@ p = 
u(x p (t),i). Hence, these values must be interpolated. In the present work, a standard fourth 



order interpolation scheme was implemented, following the recommendations of Yeung and 



Pope (1988) and Balachandar and Maxey (1989). 



In the following, the notation '@p' is adopted for arbitrary functions /(x, t), i.e., 

f @p (t) = f(x p (t),t). 



(3) 



For example uj refers to the space- and time-dependent solution of the Navier-Stokes equa- 
tions whereas Uf &p refers to the time-dependent fluid velocity seen by the particle. Corre- 
spondingly, Quf refers to the space- and time-dependent solution of the filtered Navier-Stokes 
equations whereas (£?u/) @ refers to the time-dependent filtered fluid velocity seen by the 
particle. 

Equation ([2]) is a stiff differential equation for small Stokes numbers. The numerical 
scheme for integrating equation ^ must be capable to handle this. Therefore, equation ^ 
was solved by a Rosenbrock- Wanner method (see Hairer and Wanner, 1990). This method 



is a fourth order method with adaptive time stepping. The stiff term in equation ^ is 
linearised in each time step and discretised by an implicit Runge-Kutta scheme. 

The code was validated via probability density functions (PDFs) for the particle acceler- 
ation. To this end, a DNS of forced isotropic turbulence at Re\ = 265 on 1030 3 grid points 



was conducted. This data was then compared to data from a DNS conducted by Biferale 



et al. (2004) and an experiment conducted by Ayyalasomayajula et al. (2006). Biferale et al. 



conducted a DNS at Re\ = 280 and traced inertia free particles (i.e. St — 0). Ayyalasomaya- 
jula et al. 's experiment was at Re\ = 250 with particle Stokes numbers St = 0.09 ± 0.03. 
Correspondingly, in the present simulation two particle fractions were traced, one at St — 
and another at St — 0.1. Each fraction consists of 960000 particles. Figure [2] shows that the 
results from the present simulations agree very well with the referenced data. 
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Figure 2: Probability density function -P(a) of particle acceleration a for validation of the code. X-axis is 
normalised with respect to the (Stokes number dependent) rms value of a. Triangles: reference DNS of 
St = particles conducted by Biferale et al. ( 2004|) . Squares: reference experiment of St — 0.09 ± 0.03 
particles conducted by Ayyalasomayajula et al. ( |2006 ) (renormalised). Continuous and dashed lines: present 
DNS at St = and St = 0.1, respectively. 



For model assessment, 24 fractions of particles were traced with 80000 particles per 
fraction. The particles were initialised at random positions (homogeneous distribution) inside 
the computational box and traced until a statistical steady state was obtained. Then, 1000 
time records were taken within a time span of T = 250X/u rms for computing statistics. The 
temporal resolution of the statistics equals approximately the Kolmogorov time scale. With 
this temporal resolution, the Lagrangian correlation functions could be resolved for all Stokes 
numbers. The time span was large enough to guarantee that averaging in time cancels out 
oscillations caused by the forcing scheme. 

In terms of particle time scales, T is large enough to guarantee reliable statistics. From 
e = l5vu 2 rms /\ 2 it follows that T/r p = 250y/lE/ St ~ 968.2/5*. In all simulations, St < 100, 
thus T/r p > 9.68. Hence, statistics were sampled over at least 9.68 times the particle 
relaxation time. 



4. Analysed particle-LES models 



In the present section, the three particle-LES models under consideration are presented. 
They are the approximate deconvolution method which was proposed for particle-laden flow 



independently by Kuerten (2006b) and Shotorban and Mashayek (2005) and two stochastic 
models, proposed by Shotorban and Mashayek (2006) and Simonin et al. (1993). In the 
following, the models are stated and the numerical implementation used in this work is 
explained. 



4-1. Approximate Deconvolution Method (ADM) 

Model statement. ADM is well established for incompressible single phase flows 



sec 



and Adams 


1999 


Schlatter , 


2004; 


Stolz, Adams and Kleiser 



Stolz 
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Shotorban and Mashayek (2005) and Shotorban et al. (2007) analysed the capabilities of 

ADM 
f@P 



ADM for particle-laden flow. With ADM, the fluid velocity seen by the particle u4£^ is 



computed from 

N 



nfg M = (uf M ) @p = £((X- GY Gn f ) @P = (n ADM Gn f ) @p . (4) 



n=0 



Here, X stands for identity. iV is the number of deconvolution steps. % ADM i s called 
defiltering operator because it is supposed to approximate the inverse of G- 

Equation Q is solved once per time step and the particle velocity is computed from 

J..ADM p 

QU P _ c D ne p („ADM ADM\ /e\ 

The operator H = X — G can be interpreted as extractor of subgrid scales. With this 
operator, u^ DM can be written as 

N N 

u adm = n n Gu f = H n (X - U) u f = (X - U N+l ) u f . (6) 

n=0 n=0 

For N — > oo the transfer function of 'H N+1 equals zero for the resolvable scales 
(||k|| < k c ) and one for the unresolvable scales (||k|| > n c ). This shows that for large N, the 
effect of ADM can be interpreted as improving the LES filter towards a sharp spectral filter. 

Implementation of the model in this work. In the present work, the ADM defiltering operator 
"H ADM was computed in three different ways. First, it was computed as proposed by Kuerten 



(2006b). Second, it was computed making use of the DNS spectrum and third, a model 



spectrum was used. 

If a dynamic Smagorinsky model is used as fluid-LES model, then Kuerten (2006b) 
proposes to compute % ADM as approximate inverse of the corresponding test filer. In his 



work and in the present work, this is a box filter. Kuerten (2006b) approximates its inverse 
by a second-order Taylor expansion in the filter width. The transfer function of this filter is 
shown in figure [3j In the following this approach is referred to as ADM Kuerten . 

However, it is not clear whether an inverted box filter gives highest accuracy. Therefore 
ADM was tested by two more approaches. In both approaches, the ADM filter is constructed 
such that the product of filter transfer function and LES spectrum is as close as possible 
to a target spectrum under the constraint that the filter stencil covers up to 5 3 LES cells. 
The target spectrum is either the DNS spectrum or the model spectrum proposed by |Pope| 



(2000). The results from the corresponding defiltering operators are referred to as ADM DNS 
and ADM mod , respectively. The corresponding transfer functions are also shown in figure [3j 
Evidently ADM mod leads to a very much stronger amplification around k c than ADM DNS . 
This was to be expected because around k c the model spectrum is higher than the DNS 
spectrum, cf. figure [T} 
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Figure 3: Transfer functions of the defiltering operators for the three implemented ADM approaches. 



ADM does not take explicitly into account that the model itself affects the particle path. 
More precisely, the model inherently assumes that the resolved spectrum seen by the particle 
is not modified by the model itself. In order to differentiate between this model assumption 
and other approximation errors of the model, the a priori analysis was conducted such that 
the model does not affect the particle path. 

More precisely, in the a priori analysis for ADM, for each particle two different values for 
the particle velocity were computed simultaneously. One value, referred to as DNS particle 
velocity, is the velocity obtained from the DNS flow field. The second value, referred to as 
modelled particle velocity, is the velocity obtained from a filtered DNS field and ADM. The 
particles were tracked with the DNS particle velocity and statistical samples were taken from 
the modelled velocity This approach basically tests whether ADM is capable to do what it 
is supposed to do, neglecting the effect of ADM on the particle path. 

4-2. Langevin-based models proposed by Shotorban and Mashayek and Simonin et al. 

ADM cannot reconstruct scales smaller than the LES grid. In order to circumvent this, 



Shotorban and Mashayek (2006) and Simonin et al. (1993) propose stochastic models based 



on a Langevin equation for the fluid velocity seen by a particle. Such models were originally 



developed for inertia free particles by Pope (1983), Heinz (2003) and Gicquel et al. (2002), 



referred to as generalised Langevin models. 



Statement of the model proposed by Shotorban and Mashayek (2006). Shotorban and Mashayek 



adopted generalised Langevin models for inert particles. They propose to compute the fluid 
velocity seen by the particles u^@° from the stochastic differential equation (Langevin equa- 
tion) 



du 



Sho 
f@p,i 



du 



dt 



dx,, 



dt- 



u 



Sho 
f@p,i 



dt + dWi (7) 
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and the particle velocity from 



du 



Sho 



24r n 



u 



Sho 



) dt. 



The reader is reminded that denotes 'at the particle position', cf. equation ([3]). The 
first term on the right hand side of equation ^ is the filtered material derivative of the fluid 
velocity and can be computed from the right hand side of the filtered Navier-Stokes equation. 
The second term is a drift term for the random variable Uj@ p , leading to a relaxation of u^ p 
against (Quf) mp . The last term is a diffusion term for u^ p . W denotes a Wiener process 
and e is the (modelled) dispersion of subgrid scale kinetic energy. The model parameters Tl 
and Co are specified below. 



Statement of the model proposed by Simonin et al. (1993). Simonin et al. (1993) also propose 



to model the fluid velocity seen by the particles by a stochastic process. Fede et al. (2006) 



presented in detail how to deduct Simonin et al.'s model for particle-laden flow starting from 
the Navier-Stokes equations. This results in a different Langevin equation than the equation 
proposed by Shotorban and Mashayek. 

In contrast to Shotorban and Mashayek, Simonin et al. propose to transport the resolved 
scales by particle velocity (and not by fluid velocity). The model can be formulated via a 
Langevin equation for the unresolved scales 



du 



Sim 1 
f@p,i 



(9) 



nj = Q (uiUj) — QuiQuj is the SGS stress tensor. The model constant Co is equivalent to Co 
of Shotorban and Mashayek's model. For isotropic turbulence, Fede et al. (2006) recommend 



r 



(I + Po) e 



k 



1 



i i , 



SfJS 



1 1, 



(10) 



Sij denoting the Kronecker delta function. This form of T was adopted in the present work. 



The fluid velocity seen by the particles is then computed from u 



and the particle velocity is computed from 



Sim 
f@P 



+ u 



Sim' 
f®P 



du 



Sim 



c D Rer 



2At ti 



v ( 

f@P 



(u 



u 



Sim 



) dt. 



'ID 



Implementation of the model in this work. For model closure, two parameters need to be 
specified, namely the time scale Tl and the Kolmogorov constant Co. Based on the recom- 
mendation of the model's authors, in the present work Tl was set to 



Tr 



k 



sgs 



a 



£3/2 

A 



(12) 



k sgs denotes the subgrid kinetic energy and was computed from the DNS data, e denotes the 
SGS rate of dispersion. The model constants C t and Co were set to C e = 1 and Co = 2.1, 
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following iBerrouk et all (120071) and IGicquel et all (12002b . 



For the model of Simonin et al. (1993), the SGS stress tensor r was computed in accor- 



dance with the ffuid-LES model, i.e., using an eddy viscosity hypothesis. 

As mentioned above, for the a priori analysis of ADM, the particles were traced along 
the path computed from DNS. For the stochastic models, this would be in contradiction 
to the model assumptions because the model takes explicitly into account that the particle 
path depends on the modelled small scale fluctuations (see Shotorban and Mashayek, 2005 



Fede et al. 2006). Therefore here the particle paths were computed from the modelled fluid 



velocity. 

The stochast ic differential equations (j 7|) and ^ were solved by an Euler-Maruyama 
scheme (see e.g. 



Kloeden and Platen 



2010). The stiff terms 



-uf^jT L and Ti<v 



Sim' 



were 

discretised implicitly. Shotorban and Mashayek (2005) and Fede et al. (2006) used an explicit 
Euler-Maruyama scheme. These authors focussed on small Stokes numbers. In the present 
simulations no significant differences between explicit and implicit discretisation was found 
at small Stokes numbers. At high Stokes numbers, the explicit scheme was found to produce 
significantly worse results. In particular, the kinetic energy seen by the particles explodes at 
high Stokes numbers when using an explicit scheme. It should be noted that the terms under 
consideration are linear and therefore implicit schemes do not produce any computational 
overhead. In the following, 'Sho' denotes results from the model proposed by Shotorban and 



Mashayek (2005) and 'Sim' denotes results for the model proposed by Simonin et al. (1993). 



5. A priori and a posteriori analysis of particle-LES models 

The present section contains results from a priori and a posteriori analysis for assessment 
of the three models presented in section |4| The analysis comprises the kinetic energy seen 
by the particles, particle kinetic energy and rate of dispersion. 

5.1. Assessment of ADM 

Figures [4] to [6] show the kinetic energy seen by the particles, particle kinetic energy and 
the rate of dispersion for the three ADM implementations under consideration. In addition, 
results from filtered DNS and LES without particle-LES model are shown. In order to obtain 
comparable results, the presented data from filtered DNS without model corresponds to the 
data which ADM receives, i.e., in particular in the filtered DNS the particles were traced 
along unfiltered paths. 

The kinetic energy seen by the particles (figure [4]) shows a Stokes number dependence 
due to particle clustering. The Stokes number dependence can be observed in all simulations 
although the results from filtered DNS and LES are shifted towards higher Stokes numbers 
in comparison to the unfiltered DNS. This shift is not corrected by ADM. 

As expected, the kinetic energy seen by the particles is lower in LES or filtered DNS 
than in unfiltered DNS. ADM leads to a clear improvement of the kinetic energy seen by the 
particles although a clear gap between results from ADM and DNS persists. This gap was 
to be expected because ADM does not reconstruct the smallest scales. 

Among the three ADM approaches, Kuerten's model shows least improvement. This is 
not surprising because Kuerten's approach corresponds to a single defiltering step, N = 1. 
The other two approaches must perform better here because the wider stencils allow for 
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larger values of N, leading to higher kinetic energy. On the other hand it must be mentioned 
that ADM Kuerten is computationally less expensive than ADM DNS or ADM mod and somewhat 
more generalistic in the sense that ADM Kuerten only assumes a specific filter whereas ADM DNS 
and ADM mod are based on the spectra of this specific flow. 

On first sight, the comparison of ADM DNS against ADM mod is surprising. ADM mod 
shows closer resemblance to the unfiltered result than ADM DNS although ADM D is based 
on the unfiltered field. Actually this is an effect of two errors cancelling out each other. 
Around the cutoff wavenumber the model spectrum is higher than the DNS spectrum but 
interpolation of the fluid velocity on the particle position leads to strong damping around the 
cutoff wavenumber. In other words, for ADM mod the damping properties of the interpolation 
scheme brings the spectrum seen by the particles closer to the spectrum of the DNS flow 
field. 

Figure [5] shows k u , the kinetic energy of the particles themselves. As expected from the 
observations on fc„@ p , k u is underestimated by ADM. With St — > 00 the error vanishes. 

Most interesting is the rate of dispersion, shown in figure [6| It was computed from 
the product of particle kinetic energy and integral time scale. The results from a priori 
and a posteriori analysis differ qualitatively. In the a priori analysis, ADM leads to an 
underprediction of the rate of dispersion for all Stokes numbers. The reader is reminded 
that in the a priori analysis of ADM, the particles were traced along DNS particle paths. 

The a posteriori analysis shows that actually in LES the rate of dispersion is overesti- 
mated by ADM as a result of too high integral time scales. In the a posteriori analysis of 
ADM, the particles were traced with the modelled particle velocity. Thus, the qualitative 
difference between a priori analysis and a posteriori analysis could be twofold; either due to 
the difference in particle paths or due to a defect in the fluid-LES model. Additional sim- 
ulations showed that the principal reason is that the fluid-LES model leads to too high life 
times for the large eddies and thus to an overprediction of particle dispersion. Consequently, 




St St 



— (unfiltered) DNS, --■ filtered / LES, — — ADM Kuerten , — — ADM DNS , — . — ADM mod 

Figure 4: A priori (left) and a posteriori (right) analysis of ADM, kinetic energy seen by particles. 
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it depends on the fluid-LES model whether the rate of dispersion is under- or overestimated 
by LES. 

However, a priori and a posteriori analysis show that even for high Stokes numbers where 



small scale effects should be negligible (see Yamamoto et al. , 2001), the rate of dispersion is 
not predicted correctly by ADM. 
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Figure 5: A priori (left) and a posteriori (right) analysis of ADM, particle kinetic energy. 
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5.2. Assessment of the Langevin-based models 

The testcase for the Langevin-based models is again isotropic turbulence at Re\ = 52. 
Figure [7] shows the kinetic energy of the fluid seen by the particles. In particular the results 
from the a priori analysis are very disappointing. The model of |Shotorban and Mashayek 
(2006) shows too little k u @ p for St = 0.1 although k sgs was computed from the DNS data. 

We found that this is an effect of the convective term in the model, ( Q ' 



dt 



If this term is neglected, then k u @ p is correctly predicted by the model. 

Concerning the a priori analysis, for low Stokes numbers the model of Simonin et al. 



(1993) leads to a more accurate prediction of k u © p than the model of Shotorban and Mashayek 



(2006) but for high Stokes numbers one can observe the inverse. In the a posteriori analysis, 
the model of Simonin et al. (1993) gives a satisfactory prediction for k u @ p whereas the model 



of Shotorban and Mashayek 



(2006) leads to an overestimation. 



The error in the kinetic energy of the particles is simply a consequence of the error in 
the kinetic energy seen by the particles, cf. figure |8j In the a priori analysis of the model 
of Simonin et al. (1993), the excess in k u @ p for high St is not apparent in k u . On the other 



hand, in the a posteriori analysis of the model of Shotorban and Mashayek (2005), the excess 
in k u @ p for high St is also visible in k u . This is probably due to the different time stepping. 
In LES, the time step size is higher than in DNS. For the particles St = 100, the ratio 
between particle relaxation time and time step size is about 1000 for DNS and about 200 
for LES. Thus, the St = 100-particles can rather follow the modelled fluctuations in LES 
than in DNS. Therefore, the excess in k u @ p at high Stokes number is rather reflected in the 
a posteriori analysis than in the a priori analysis. 

Concerning the rate of dispersion, figure [9j the result from the a priori analysis is very 
discouraging, the a posteriori results are somewhat better concerning accuracy of the models. 
Nevertheless, these tests show that both models do not necessarily improve the result of the 
LES in comparison to an LES without particle-LES model or LES with ADM. In particular 
at high Stokes numbers it might be recommendable to use no model instead of one of the 




Figure 7: A priori (left) and a posteriori (right) analysis of the Langevin-based models, kinetic energy seen 
by particles. -, q 



stochastic models, in accordance with the numerical results of Shotorban and Mashayek 
([2005]). 
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6. Conclusions 

We have presented data from DNS, filtered DNS and LES of particle-laden isotropic 
turbulence with three different models for the effect of the unresolved scales on the parti- 
cles (particle-LES models). The models under consideration are approximate deconvolution 



(ADM) as proposed by Kuerten (2006b) and two stochastic models, based on the works of 



0.8- 



0.6 h 
0.4- 
0.2- 



?0" 



a priori analysis 




- unfiltered 

- filtered 

- Sho 

- Sim 

jj i 



10° 10 1 
St 



10 2 



Q. 

XL 




Figure 8: A priori (left) and a posteriori (right) analysis of the Langevin-based models, particle kinetic 
energy. 
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Shotorban and Mashayek (2006) and Simonin et al. (1993). The present work allows for the 



first time a direct comparison of these models because all models were assessed on the same 
testcase. 

The models were assessed by a priori and a posteriori analysis. The analyses comprise 
the kinetic energy seen by the particles, particle kinetic energy and rate of dispersion. 

ADM was implemented in three different ways, one following Kuerten (2006b) and two 
approaches via optimisation against spectra. With all three approaches a defect in the rate 
of dispersion could be observed. ADM was found to underestimate the rate of dispersion if 
applied on the filtered DNS field. In LES however, ADM was found to overpredict the rate 
of dispersion. This discrepancy was explained as an effect of approximation errors in the 
fluid-LES model. 

The stochastic models showed very poor performance in the numerical simulations. In 
particular for high Stokes numbers, LES or filtered DNS with stochastic models showed 
larger difference to DNS results than LES or filtered DNS without particle-LES model. 

According to these results, the stochastic models are not recommendable because, in de- 
pendence of the configuration, LES without particle-LES model can perform better, some- 
times even tremendously better, than LES with a stochastic particle-LES model. On the 
other hand, results from ADM are quite promising. At least ADM was found to lead to an 
improvement for all Stokes numbers. 

On the other hand, ADM only enhances the resolved scales but does not actually model 
scales which cannot be represented on the grid. At high Reynolds numbers, where the LES 
grid is very coarse due to computational limitations, ADM can be expected to perform worse 
than at low Reynolds numbers. Therefore a new particle-LES model is needed for LES at 
high Reynolds numbers. 
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